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Abstract. We study the Bloch oscillations (BOs) of two-component Bose-Einstein 
condensates (BECs) trapped in spin-dependent optical lattices. Based on the derived 
equations of motion of the wave packet in the basis of localized wave functions of the 
lattice sites, the damping effect induced by the intercomponent and intracomponcnt 
interactions to the BOs is explored analytically and numerically. We also show that 
such damping of the BOs can be suppressed entirely if all the atom-atom interactions 
are modulated synchronously and harmonically in time with suitable frequency via the 
Feshbach resonance. When the intercomponent and the intracomponent interactions 
have inverse signs, we find that the long-living BOs and even the revival of the BOs can 
be achieved via only statically modulating the configuration of optical lattices. The 
results provide a valuable guidance for achieving long-living BOs in the two-component 
BEC system by the Feshbach resonances and manipulating the configuration of the 
optical lattices. 
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1. Introduction 

The system of ultracold atomic gases in optical lattices has become a nice experimental 
platform to simulate the effects in condensate matter physics pp. The well controllability 
of such system makes many sophisticated effects in condensed matter physics be well 
studied in this system [2]. The Bloch oscillation (BO) is the oscillatory motion of a 
quantum particle in a periodic potential when it is subjected to an external force. It 
was originally predicted in solid-state system where the motions of the electrons in 
tilted periodic potentials undergo coherent oscillations [3]. The formal resemblance 
between electrons in crystals and Bose-Einstein condensates (BECs) in optical lattices 
has inspired an extensive interest to explore the BOs in optical lattice system. The 
successful observations of BOs have been reported for atoms in interacting BECs 
[U El El [7] and for ensembles of noninteracting quantum-degenerate fermions [8] in 
tilted optical lattices. 

However, the perfect BOs can only be available in the ideal case where there is 
no interactions among the atoms of BECs. Practically, due to the intrinsically weak 
interactions of atoms, the momentum distribution of the BECs will show a rapid 
broadening, which causes the atoms of BECs to lose their phase coherence, i.e. the 
dephasing effect [91 [10] . Consequently, the BOs in BECs cannot persist on for a long 
time. In the framework of the mean-field treatment, the motion of the BECs can be well 
described by the so-called Gross-Pitaevskii equation (GPE) with a nonlinear term. It is 
believed that such nonlinearity induced by the atom-atom interactions generally leads 
to a breakdown of the BOs, as recently studied experimentally [UJ and theoretically 
[T2~l ITB"] [H]. Therefore, a natural question is: Is it possible to prolong or even stabilize 
the BOs by some active control ways? 

Addressing this question, some progress, especially in single- component BECs, has 
been made. It has been found that a long-living BO can be induced by properly designing 
the spatial dependence of the scattering length (15] and the configuration of the optical 
lattices [16J. Gustavsson et al. [H] showed the experimental evidence that the dephasing 
time of the BOs can be much enhanced by decreasing the interaction strengths via the 
Feshbach resonance [17]. However, in many situations, such finite enhancement to the 
dephasing time of BOs is not enough, and one is desired to preserve the BOs forever. 
Recently, Gaul et al. have reported that a persistent BO of single- component BECs can 
be obtained by modulating interaction harmonically in time with suitable frequency and 
phase [TO], [18] , which can be easily realized by means of Feshbach resonance. 

So far, most of the studies of the stabilization control to BO of BECs in optical 
lattice are based on the single-component BEC case. Compared with single-component 
BEC, the two-component BEC system may exhibit more novel physical effects due to 
the condensate mixtures [191 [2^ EH [22]. In this system, both intracomponent and 
intercomponent atom-atom interactions take effects upon the nonlinear behavior of the 
BECs. The two-component BECs can be experimentally realized by spin-dependent 
optical lattices, which hold the BECs with the two components composed of two distinct 
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hyperfine states of the same atomic species [23] . 

In the present paper, we study systematically the modulation stabilization of the 
BOs of two-component BECs in optical lattices. We mainly use two control ways, one 
is by modulating periodically the interactions via Feshbach resonance; the other is by 
tuning statically the parameters of the optical lattices. We will show that the stable BOs 
can be obtained when the interactions are modulated synchronously and harmonically in 
time with suitable frequencies. Moreover, if the intercomponent and the intracomponent 
interactions have inverse signs, the long-living BOs and even the revival of the BOs can 
be achieved via only tuning the the relative separation between lattices. 

The paper is organized as follows. In Sec. II, we discuss the methods and formalism 
used in this work. In Sec. Ill, we explore quantum manipulation of BOs from two aspects 
according the signs of the intracomponent atom-atom interactions. Finally, a summary 
is given in Sec. IV. 

2. Model and formulation 

2.1. Gross-Pitaevskii equations for the two- component BECs in an optical lattice 

We consider two-component BECs which are composed of bosonic atoms of the same 
isotope but having different internal spin states, e.g. 87 Rb atoms in hyperfine states 
\F — 2, niF — 2) and \F — l,mF — — 1) [21] • The BECs are trapped in spin-dependent 
optical lattices. The dynamics of the system is governed by the coupled GPEs under 
the mean field approximation, 

= [ - ^v 2 + Vi+ E^(*)l*il a ]*<> (!) 

where $j (i = 1, 2) is the macroscopic condensate order parameter of the i-th component 
with identical mass M. The time-dependent interaction coefficients are given by 
gij(t) = 4irh 2 aij(t) /M with Oy(t) being the s-wave scattering length which can be 
controlled via the Feshbach resonance induced by the modulated magnetic field . The 
external potential felt by the i-th component can be decomposed into Vi = V c +VLi, where 
V c = fz is a linear potential induced by a constant force / and Vu is trapping potential 
of the optical lattice. The additional weak potential V c tilts the optical potentials and 
drives coherent oscillations (5] . The trapping potentials for different components can be 
explicitly written as Vl\ = U p sin 2 (k L z + |) and V^ 2 = U p sm 2 (kLZ — |), where U p is the 
depth of the ID lattice potentials, ki is the wave vector of the lasers used to construct the 
optical lattice, and 9 is the polarization angle of the two counterpropagating laser beams 
to form the standing wave configuration of the optical lattice (2Ql [23]. By changing 9, 
one can also control the separation between the two potentials. 

When the linear field is too weak to induce Landau-Zener tunneling [25l [26], BO 
can be described by an adiabatic evolutions of the BECs in the lowest lattice band. In 
collective coordinates [9], the condensate order parameter $i(r, t) can be expanded as a 
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linear combination of the wave packets localized at the individual lattice sites, i.e. the 
Wannier wave functions d) n .[r), as 



Hr,t) = y/Ni y £ l iI>i,rH(t)K(r), (2) 

m 

where iVj is the total number of particles of the z-th component and the Wannier 
wave function (p ni satisfies the orthogonality condition / dr(fi ni (j) ni ±i = 0, and the 
normalization condition / drcf^. = 1. ^ jTli = yj iV^ (t) / NiC ldl ' n i ® , where Ni m (t) and 
^i,nj (£) are the number of particles and phase, respectively, is the amplitude of the i-th 
component trapped in the site n*. In the following, we assume that the two components 
have the same total number of particles, i.e., Ni = N 2 . Substituting Eq. (j2J) into Eqs. 
(II]), we can discretize Eqs. (JTJ into a set of coupled nonlinear equations with respect to 
different lattice site rij, 

• ; Yi,Tli — l "t" Yi,m+1 . r a / ,\ i , |2 

i^i,rH = 2 + + A "(*)l^wl 

,ni+T—l \ 2 Mi,ni, (3) 

where the overdot denotes the time derivative and i ^ j labels the two different 
components of BECs. e ijUi = jj j dr[ {hV ^ + V^.] with J = - J rfr[^V0 ni V0„ l+ i + 
0n,^0n,+i] being the tunnel parameter. = g "^J JV ' / G?r<^. describes the 

intracomponent atom-atom interaction strength. Ay = 3 ' J *J — multiplying with 
r] T = J dr(j)n.<f>n-+ T an d 'Tr-i — / ^ r 2 4 2 i+r _i describe the intercomponent atom-atom 
interaction strengths, where r] T and r/ T _i stem from the overlaps of the wave functions 
of the two components, r (in the lattice unit) is determined by the relative separation 
between the two nearest neighboring spin-dependent potentials, which can be controlled 
by the polarization angle 9. It is noted that not only all the interactions are time- 
dependent but also intercomponent nonlinear interactions depend on r] T and f] T -i- In 
Eq. (J3J) the time has been rescaled to be dimensionless as t —> ht/2J. 

In ID optical potentials, we can denote the Wannier wave function as <fi ni (r) = 
(f>(x,y)(f> ni (z). The transverse wave function can be expressed as a 2D Gaussian profile 
4>{x,y) = (po(x)(po(y), where (f>o(a) = ^ , — exp(— ■^- T ) with a a being the Gaussian 
widths in the a = x,y directions. The wave function along the direction of the optical 
lattice can be denoted as <fi ni (z) = <p (z — riid) with d = n/k L being the lattice constant. 
Based on the variational ansatz for <f>m(z), a minimum energy can be obtained when the 
width of 4>m( z ) equals to a z = ~dU p / E rec , where E rec is the recoil energy [25]. Under 
this consideration, the parameters in Eq. ([3]) can be determined explicitly as 

e 2 u p e 2 u p 

ei, ni = uni H — , e 2 ,„ 2 = um 2 — L , 



. . , 1 gu(t)Ni , r 2 d 2 , 

^ ] ^ (2^ ' ^ = exp(-— ), 

J= exp(4){^^i-^[l + exp(^i)]}, (4) 
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where u = |j and / corresponds to the weak atomic gravity. From Eq. ([3]) and the 
canonical equation ipi = g^^y, the Hamiltonian functions can be obtained 

+ ^|^,nJ 2 + A 4J (t)(r/ T |^,„ l+T | 2 

+ 1 ^,n i+ r-l H] \*Pi,nf}, (5) 

both the Hamiltonian and the norm iV^nJ 2 = are conserved. 
2.2. The dynamics of wave packet: Bloch oscillations 

In order to analyze how the interactions affect the BOs of the two-component BECs, 
we parameterize the Gaussian profile wave packet for i-th component as [H] 

li 2 



^i,m = V^expf h ipi(rii - &) + ^(n, - &) ]. (6) 



where Ki = is a normalization factor. The Gaussian wave packet is parameterized 
by four types of quantities: the center-of-mass position £i(t), the width of the wave 
packet described by ji(t), the linear phase pi(t) describing the group velocity of the 
wave packet, and the quadratic phase 5i(t) over the wave packet. The latter phase 
allows us, on the one hand, to describe the linear evolution of the wave packet for which 
the quadratic dispersion in momentum space directly translates into a quadratic phase in 
real space. On the other hand, the nonlinearity due to the atom-atom interactions also 
leads to a quadratic phase since the density near the Gaussian maximum is quadratic 
[I]. Such Gaussian profile wave packet was used to explain successfully the BO in 
Anderson-Kasevich experiment [5]. 

The dynamical evolution of the wave packet can be obtained by a variational 
principle from the Lagrangian £j = i^Z ni ipi, ni '4 , in- ~ After some algebra, the 
Lagrangian can be achieved 

r i ^ j_ -w A *(*) 

Q = PiKi 3- + e Xl cos ^ - — -= Vi 

- A^-^e"^ + r/.-ie-^- 1 ], (7) 

\/7T 



where X i = 5^ + = K t J dn^ exp[ 2(n 'J i)2 ], k = ^ with 7 2 = 7l 2 + 7 |, and 

/i r = k 2 £ 2 with £ r = & — £j + t. It is noted that in our calculation the summation 
over rij has been replaced by integration when the widths 7, are not too small [9j [22] . 
The equations of motion of the collective coordinates can be obtained from the Euler- 
Lagrange equations 

ii = e" Xl sin ft, (8) 
ji = lAe~ Xl cos pi, (9) 
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Pi = 2 * 3 y V ere"^ + •n r . 1 i T . X e-^\ ~ CO, (10) 



i = ( l_ i | )e -» cosP( + !^W + ^!^MW, 

where z = r] T (j 2 - 4£^)e'~' iT + r/ T „i(7 2 - 4£^_ 1 )e~' tT - 1 . 

To highlight the essential physics, from the coupled Eqs. fl8|)-f lTT]) . the equation of 
motion of the center of the wave packet can be recast into 



where 



£ + = (12) 

(3i = Pie~ Xi cos (14) 

It is noted that Eq. ( IT2|) can recover the equation of motion of the wave-packet center 
for single component BEC under the condition Ay = as [9] 

6 + o,6 + co% = ^[7^(0) - 7^)], (15) 

which is a standard equation of motion for a harmonic oscillation with an effective 
damping rate a(t). Under the ideal condition An = 0, a(t) = and the system 
undergoes a perfect oscillation, which is the perfect BOs with the frequency w as a 
result of the driven field. The damping of the BOs is caused by the intracomponent 
interactions An, which means that the nonlinearities inevitably lead to the breakdowns 
of the BOs. In the situation of two-component BECs, the intercomponent interaction 
Ay will add another nonlinearity to each individual component of the BECs, which 
gives an additional contribution to the damping rate, as shown in Eq. f|T3|) . Besides 
the damping rate, the intercomponent interaction also exerts an effective driven force 
to BOs. Explicitly, under the condition that Ay is small compared with the constant 
force / felt by the BECs, one can rewritten Eq. ( 1T21) as 

ii + ctiii + u 2 £,i = r]i{Aii, Ay), (16) 

which shows that the perfect BOs of the system are distorted by the cooperative action of 
the effective damping rate a(t) and the effective driven force r}i{An, Ay). It is noted that 
the frequency u of the BOs, which corresponds to the inverse of the right-hand side of 
Eq. (JTUjl . for the two-component case is slightly detuned from uj by the intercomponent 
interactions Ay. The explicit analysis of the BOs governed by a(t) and rji(An, Ay) will be 
shown by the quantitative calculations with the experimentally adjustable parameters 
in the next section. 



3. Quantum manipulation of Bloch oscillation in the optical lattice 



From the above analytical results, we can see that the nonlinearities contributed from 
both of the intercomponent and intracomponent interactions generally lead to the 
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breakdowns of the BOs of the BECs in the optical lattices. Therefore, the interaction 
coefficients of such nonlinearities, which are essentially determined by the s-wave 
scattering lengths ay, have profound impact on the dynamics of two-component BECs 
in the optical lattices. The magnitude and sign of these parameters sensitively influence 
dynamical behaviors of this ultracold boson system. In cold-atom experiments, Feshbach 
resonance is a quite effective mechanism that can be used to modulate g,ij. Inspired 
by a recent experimental investigation of the ultracold molecule production via a 
sinusoidal magnetic field modulation to the interaction coefficient around the Feshbach 
resonance [27], we intend to explore the possibility to stabilize the BOs via such 
periodic modulations to the interaction coefficients g^ in the following. Besides the 
magnetic field induced Feshbach resonance, another way to modulate the dynamics of 
the BECs in our system is via adjusting the configuration of the optical lattices. The 
separation between the spin-dependent potentials felt by the two components of BECs, 
which can be adjusted by the polarization angle 9 of the lasers, essentially determines 
the intercomponent interactions of the BECs. We also examine the influence of the 
separation on the dynamics of BOs. The combined effect of the Feshbach resonance and 
a periodic external potential has been widely studied [2EJ EHJ 130] . 

3.1. a u (0) > 0,oy(0) > 

In this case, both of the intercomponent and intracomponent interactions are repulsive. 
Without loss of generality, we assume that the intercomponent interactions Ay relate 
to the intracomponent interactions as: A 12 = A 2 i = \M-n-A-22 i n our numerical 
simulations. 

To see the effect of nonlinearities on the BOs, we plot in Fig. [1] the attenuations 
of BOs without modulations. From the time-dependent behaviors of the wave-packet 
centers & [Fig. GJa)] we can see clearly the breakdown of BOs with time evolution. 
Such damping oscillations are manifested by the behaviors of the effective damping 
rates a« [Fig. H^c)] and the driven force rji [Fig. CDT)], where ccj are always positive 
and increase with time, and r\i tend to constants after the same cycles as Compared 
with the case for the single-component BECs [9], the presence of the intercomponent 
interactions, sharing the same sign with intracomponent ones here, play the role as an 
additional nonlinearities and speed up the collapses in our two-component situation. So, 
the stronger the intercomponent interactions are, the faster the BOs damp. Fig. QJd) 
shows that Pi, the associated momenta of £j, increase linearly with time. This can be 
understood from the analysis of Eq. ( flOl) . The first term of the right hand of Eq. (llOjl . 
which is contributed from the intercomponent interactions, is much smaller than the 
second term, which is contributed from the linear potential. Consequently, the time- 
dependent behaviors of pi are dominated by the second term, i.e. Pi{t) ~ Pi(0) — ut. 
Fig. QXb) shows that the widths 7$ undergo breathing oscillations and soon approach 
constants. The associated momenta 5i are also divergent, as shown in Fig. [He). All these 
time-dependent behaviors indicate that the system is set into a macroscopical quantum 
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Figure 1. (Color online) Attenuations of BOs without modulating the interactions. 
Solid and dashed lines indicate two individual components with the intracomponent 
interactions as An = 20 and A22 = 15, respectively. The other parameters are 
chosen as: U p = 16E rec , ui = 2, and r = 0.5. The initial conditions are set as: 
pi(0) = P2 (0) = 0, <5i(0) = S 2 (0) = 0, 6(0) = 1, 6(0) = 0, 7l (0) = 72(0) - 10. 




Figure 2. (Color online) The stabilization of BOs by modulating the interactions. All 
interactions are modulated by cos(wt). The same parameters and notation are used as 

Fig. m 
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self-trapping mode [3TJ E2J E3] due to the nonlinearities from the intercomponent and 
intracomponent interactions. 

Now we focus on how to stabilize the BOs in our system. We mainly use the way 
by modulating periodically the interactions via a magnetic field. 

In Fig. [2] we plot the time-dependent behaviors of the wave-packet variables under 
a cos cut, where u is the frequency of BOs, modulation to all the interactions. We 
can see that the damping of the BOs can be fully stabilized by such modulation and 
perfect oscillations is thus obtained. Such persistent phenomena can be also explained 
straightforwardly by studying the time- dependence of the damping coefficients oti [Fig. 
EJc)] and the effective driven forces rji [Fig. 12(f)]- In contrast to the positivity in the full 
evolution range in Fig. [1(c), aj in this case exhibit periodic oscillations between positive 
and negative values with definite amplitudes, which characterizes well the lossless BOs 
of £j. It means that the effective damping coefficients with alternate signs inspirit the 
system itself to guarantee the stabilizations of the BOs. Besides, we find that the driven 
forces in Fig. |2(f) are replaced by constant values completely after such modulation. 
It has been proven that there is a family of stable solutions in terms of collective 
coordinates in single- component BEC system when all the interactions are modulated 
by cos(cjt) [W\. In fact, in two-component ones, the coupled terms in the Eqs. (OQ) can 
be regarded as additional nonlinearities, which possess the same time dependence as 
intracomponent interactions so long as modulating all interactions harmonically in time 
with the same suitable frequencies. As viewed from the mean field, the two components 
can be reduced into two independent single ones. So it is understandable that such 
modulation stabilization also presents in our two-component BECs system. 

It is noted that the effect of the modulation to the BOs is sensitively dependent 
of the forms of modulating field we used. In Fig. [31 we plot the numerical simulation 
when all the interactions are modulated by a sin(a;t) field. It is found that the damping 
of the BOs manifested by & is not suppressed and the BOs is destroyed in several 
rounds of oscillation. The wave-packet widths 7, reduce their amplitudes quickly. So 
the modulation takes no effect in this case. 

In the present case, the intercomponent interactions share the same sign as 
the intracomponent ones, so the tuning of the relative separation r does not have 
constructive action to suppress the damping of BOs. However, the things are 
changed dramatically when the intercomponent interactions have opposite sign to 
intracomponent ones, as discussed in the following. 

3.2. Oii(0) > 0,^(0) < 

In this case, the original intracomponent interactions are repulsive, while the original 
intracomponent ones are attractive. Different to the above case, there are many 
interesting effects induced from the competition between such two kinds of atom- 
atom interactions. For example, the stability of static solitonic excitations in two- 
component BECs have been analyzed within the Gross-Pitaevskii approximation [34]. 
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Figure 3. (Color online) The attenuations of the BOs with modulating interactions 
by sin(wt). r — 0.1 and the same other parameters and notation as Fig. [T] 

For convenience, we assume the intercomponent atom-atom interactions to be A 12 = 
A 2 i = — \M-n-A-22 i n our numerical simulations. 

As analyzed in above case, the BOs are destroyed by the nonlinearities. From this 
point, the dynamical behaviors of the wave packet in the present case show no difference 
to the above one. However, we can prolong the coherent time of the BOs dramatically 
by tuning the relative separation r of the potentials felt by the two components in the 
present case. To confirm this, we plot the time evolutions of the wave-packet centers 
£i and widths 7$ in Fig. H] for different relative separations. A large r means a large 
distances between the nearest neighbors of the Wannier wave functions, which in turn 
induces a small intercomponent interaction rate t] T . Fig. H^a,d) show the breakdown of 
the BOs when the intercomponent interactions are small for a large r. With the decrease 
of r, the intercomponent interactions get stronger. The damping of BOs are obviously 
slowed down [Fig. I3|b,e)]. Especially, it is noticed that 7$ shows revival at about t = 150, 
as the dashed line in Fig. S(e). It provides a valuable guidance that the dynamics 
of the system would show revival in this case. If the relative separation r is further 
reduced so that the two lattices are extremely close, the BOs show obvious revival [Fig. 
Hfc,f)]. This phenomenon is caused by the competition between the intercomponent and 
intracomponent interactions. Because the intracomponent interactions have inverse sign 
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Figure 4. (Color online) Attenuations of BOs without modulating interactions for 
different potential separations r. Solid and dashed lines indicate two individual 
components with the intracomponent interactions as An = 20 and A22 = 18, 
respectively. The other parameters are chosen as: U p = l6E rec , lu — 2, and r = 0.5 
for (a, d), r = 0.1 for (b, e), and r = 0.05 for (c, f). The initial conditions are set as: 
Pl (0)=p 2 (0)=0, 5 X (0) = S 2 (0) = 0, £i(0) = 6(0) =0, 7l (o) =72 (o) = 10. 



with the intercomponent ones, the nonlinearities contributed from the intracomponent 
interactions are partially counteracted by the intercomponent ones. For the full overlap 
at r = 0, the two components perfectly mix together and attenuations of BOs reappear. 
To sum up, the coherent time of the BOs can be much enhanced by only tuning the 
relative separation r of the optical lattices. 

However, in many situations, such finite enhancement to the dephasing time of BOs 
is not enough, and one is desired to preserve the BOs forever. This actually can also be 
achieved by modulating the interactions by a suitable time- dependent field. Fig. |5]shows 
that the BOs are entirely stabilized by modulating interactions harmonically in time 
with the same frequency as the one of the BOs. Such stabilization is independent of the 
magnitudes of the nonlinearities, so the behaviors for different r under the modulation 
are same, as shown in Fig. |5]^a,b). 

4. Conclusions 

In summary, we have studied analytically and numerically the dynamical behaviors of 
the BOs for two-component BECs trapped in combined potentials consisting of linear 
potentials and spin- dependent optical lattices. We found that the damped BOs can 
be stabilized when all the atom-atom interactions are modulated synchronously and 
harmonically in time with Bloch frequency. Moreover, if the intercomponent and the 
intracomponent interactions have inverse signs, it has been shown that the dephasing 
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Figure 5. (Color online) The stabilization of BOs for different relative separations r 
under the cos(cjt) modulation, r = 0.1 for (a, c) and t = 0.5 (b, d). The same other 
parameters and notation as Fig. [4] 

time of BOs can be much enhanced by decreasing the relative separation of the two 
potentials felt by the two components. Our results provide a valuable guidance for 
achieving long-lived BOs in the two-component BEC system by the Feshbach resonances 
and manipulating the configuration of the optical lattices. 
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